Dynamical Effects and Phase Separation in Thin Films 
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We study phase separation in thin films using the Navier-Stokes Cahn-Hilliard equations in the 
lubrication approximation, modeling substrate- film interactions with a van der Waals potential. We 
investigate the thin-film equations numerically and compare them with experimental results. We 
find that the model captures the qualitative features of real phase- separating fluids, in particular the 
tendency of concentration gradients to produce film thinning and surface roughening. The ultimate 
outcome of the phase separation depends strongly on the dynamical backreaction of concentration 
gradients on the flow, as we demonstrate when a shear stress is applied at the film's surface. When 
the backreaction is small, the phase domain boundaries align with the direction of the imposed 
stress, while as the backreaction is made larger, the domains begin to align in the perpendicular 
direction. 
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When a binary fluid is cooled below the critical tem- 
perature, the homogeneous state is energetically un- 
favourable and the system spontaneously phase-separates 
and forms domains rich in either fluid component p], Q • 
Due to the relevance of phase-separating thin films in in- 
dustrial applications [3j], many experiments and numeri- 
cal simulations focus on understanding how phase sepa- 
ration is altered if the binary fluid is confined in a thin 
layer. We propose a lubrication approximation based 
on the coupled Navier-Stokes Cahn-Hilliard equations 
to explain the main features of these studies. 

Several recent experiments have clarified the different 
regimes of domain growth in a binary thin film. Wang 
and Composto [4| have identified early, intermediate, and 
late stages of evolution. The early stage comprises three- 
dimensional domain growth, while the intermediate stage 
is characterized by the formation of wetting layers at the 
film boundaries, the thinning of the middle layer, and 
significant surface roughening. Due to the thinning of 
the middle layer, the sandwich-like structure breaks up 
and matter from the wetting layer flows back into the 
bulk. Thus, a late stage is reached, consisting of bubbles 
coated by thin wetting layers. This characterization of 
the evolution has been seen in other experiments [H, Q , 
although clearly a variety of behaviors is possible, de- 
pending on the wetting properties of tnhe mixture. Our 
model captures the essential features of this evolution, 
in particular the tendency for concentration gradients to 
promote film rupture and surface roughening. 

In a series of papers, Das et al. @, 0] investigate the 
behaviour of binary fluids with wetting. In [TJ they spe- 
cialize to ultra-thin films. In bulk mixtures, where one 
of the fluid components is preferentially attracted to the 
boundary, a layer rich in that component may be estab- 
lished there, followed by depletion layer, and so on. This 
so-called spinodal wave propagates into the bulk Q. In 



ultra-thin films, the film thickness is less than a single 
spinodal wavelength and the spinodal wave is suppressed. 
Two distinct outcomes of phase separation are identified, 
depending on whether one binary fluid component wets 
the film boundary completely or partially. Our focus will 
be on the partially wet case. In this wetting regime, both 
fluid components are in contact with the film boundaries. 
The authors find an ultimate state of domain formation 
extending in the lateral directions and growing in time as 
t 1 / 3 , a result that indicates domain growth by Lifshitz- 
Slyozov diffusion Q. 

These papers elucidate the role of wetting and film 
thickness on the process of phase separation, although 
they do not discuss hydrodynamics or the effect of free- 
surface variations on domain formation. In this paper, 
we therefore focus on ultra-thin films with a variable free 
surface, and for simplicity we restrict our attention to 
the case where both fluids experience the same interac- 
tion with the substrate and free surface. The model we 
introduce is based on the Navier-Stokes Cahn-Hilliard 
(NSCH) equations [10[ and gives a qualitative explana- 
tion of these studies, in particular the tendency of domain 
formation to cause film rupture and surface roughening. 
With an applied external forcing, the model illustrates 
the salient effect of the dynamical backreaction of con- 
centration gradients on the flow, a useful result in appli- 
cations where control of phase separation is required [ll| . 

In full generality, the equations we study are 
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is the stress tensor, p is the fluid pressure, <j) is the body 
force potential and p is the constant density. Addition- 
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ally, v is the kinematic viscosity, (3 is the mixture free 
energy per unit mass, D is the Cahn-Hilliard diffusion 
coefficient, and y/j is the thickness of domain bound- 
aries. The concentration boundary condition for Eq. (pQ) 



is n • Vc = n 
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0, where n is a vec- 



tor normal to the boundary, while the velocity bound- 
ary conditions on the velocity and stress tensor are stan- 
dard [12]. We nondimensionalize these equations by us- 
ing the vertical length scale /io, the horizontal or lateral 
length scale A, and the diffusion time X 2 /D. If the pa- 
rameter e = ho/X is small, a lubrication approximation is 
possible [IH • We take the following dimensionless groups 
to be of order unity, 
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where Re is the Reynolds number, C n is the Cahn num- 
ber [TO] which provides a dimensionless measure of do- 
main wall thickness, r is a dimensionless measure of the 
backreaction strength, and C _1 is a dimensionless mea- 
sure of surface tension corresponding to the dimensional 
surface tension cr . Using these scalings, we expand the 
nondimensional version of Eq . fl]) in powers of e, follow- 
ing the method outlined in [12[, and obtain equations 
for the free surface height h(x,y,t) and concentration 
c(x,y,t), 
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where 



v = ~^ v i^ + ( x > y> h ( x > y> *)) + r ( v ^ c ) 2 > 

/i = c 3 -c-C^V ± -(W ± c). 

Here = (d x ,d y ) is the gradient operator in the lat- 
eral directions, a is the surface tension, <\> is the body 
force potential, u is a vertically-averaged velocity, p as a 
vertically-averaged pressure, and \i as the chemical po- 
tential. While the equations do not allow for vertical 
variations in concentration, we show in what follows that 
the model reproduces the qualitative features observed in 
thin binary fluids, especially in the case where both bi- 
nary fluid components interact identically with the sub- 
strate and free surface [7|. 

For thin films with h = 100-1000 nm 0, @|, the 
dominant contribution to the potential is due to van 
der Waals interactions [H, Qjl , and following these au- 
thors we take <p = A/h 3 , where A is the dimensionless 
Hamaker coefficient. To prevent rupture [12[, we study 
films where A < 0, and take A to be independent of 



the concentration level, so that both binary fluid com- 
ponents are attracted equally to the substrate and free 
surface boundaries. In this case, Eq. ([2j) possesses sim- 
ple one-dimensional equilibrium solutions, obtained by 
setting u = V_L/i = 0. From Fig. [T] we see that the one- 
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FIG. 1: (Color online) Equilibrium solutions of Eq. J2} for 
C = Cl = \A\ = 1 and r = 0.1, 1, 10,50. In (a) the valley 
deepens with increasing r although the film never ruptures, 
while in (b) the front steepens with increasing r. 

dimensional equilibrium solution of Eq. ([2]), with bound- 
ary conditions h (±00) = 1, c(±oo) = ±1, consists of a 
step-like profile for the concentration, corresponding to a 
pair of domains separated by a smooth transition region. 
Across this transition region, the height field dips into a 
valley. While the valley increases in depth for large back- 
reaction strength r, the film never ruptures. This result 
follows from the inequality h" (0) > 0, since x = is 
a local minimum. Thus, from the equilibrium condition 
u = 0, 



< 



ur' <or 
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In this way, the repulsive van der Waals potential has a 
regularizing effect on the solutions. 

Physically, the formation of the valley arises from 
the balance between the van der Waals and backreac- 
tion effects. From the solution in Fig. [TJ the capillary 
force F cap = — rd x (d x c) 2 and the van der Waals force 
F vc iw = \A\ d x h~ 3 always have opposite sign. The repul- 
sive van der Waals force acts as a non-linear diffusion [lij 
and inhibits rupture, and therefore F cap promotes rup- 
ture, a result seen in experiments [J]. The valley in the 
height field represents a balance between the smoothen- 
ing and the rupture-inducing effects. 

As in ordinary Cahn-Hilliard dynamics [2], the one- 
dimensional equilibrium solution hints at the late-time 
configuration in higher dimensions. Thus, we expect 
the multidimensional solution to comprise concentra- 
tion domains with a height field of peaks and valleys, 
with valleys occurring at domain boundaries. We have 
verified with numerical simulations that this is indeed 
the case. By using a measure of domain size (L x ,L y ) 
based on the Fourier transform of the correlation func- 
tion (c (a;, t) c (x + r, t)) [15[, we have found that the do- 
mains grow in time as t 1 / 3 , the usual Lifshitz-Slyozov 
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growth law @. Here x = (x,y) denotes the lateral coor- 
dinates and (...) denotes the spatial average. The modi- 
fied growth exponent due to hydrodynamic effects @, [l6[ 
is not observed, a result that emerges from the non-linear 
diffusive character of the height equation, which damps 
any undulations not caused by concentration gradients. 
The surface roughness arising from the concentration gra- 
dients is similar to that observed in the one-dimensional 
case and has been seen in several experiments [3, 03 . 

The dramatic effect of the reaction of the concen- 
tration gradients on the phase separation is apparent 
when we apply a surface tension gradient across the film. 
Physically, this can be realized by differential heating of 
the surface [Hj], although a surfactant will also induce 
stresses at the surface [19]. We set a = aosin/cx, where 
k = (27r/L)m = kom is the spatial scale of the surface 
tension variation and m is an integer. Then the velocity 
that drives the system becomes 



u = \h (kao cos kx, 0) 
+ \h 2 V A 



(3) 



This velocity field may also be obtained by imposing a 
shear stress r at the surface, provided r = Va [20]. We 
carry out simulations using Eq. ([3]) on a 128 x 128 grid. 
The results do not change upon increasing the resolution. 
We choose C n so that domain boundaries are resolved. 
The other parameter values are indicated in the caption 
to Fig. H 

This choice of velocity field leads to control of phase 
separation in the following manner. For small values of 
the backreaction strength, with r — > 0, the height field 




FIG. 2: (Color online) The concentration field for C = —A = 
1. Across the first row, r = and (a) t = 500; (b) t = 3750; 
(c) t = 7500; (d) t = 30000. Across the second row, r = \ 
and (e) t = 500; (f) t = 3750; (g) t = 7500; (h) t = 30000. 
The surface tension gradient is parallel to the arrow and a = 
cr sin(/cx), (Jo = 20 and k = 4/c . In Figs, (a)-(d) with r = 
0, the domains align along the arrow, while in Figs, (e)-(h) 
with moderate backreaction strength, the domains align in a 
direction perpendicular to the arrow. 

quickly aligns with the surface tension profile as in Fig.[3j 
since the strong effect of the van der Waals diffusion de- 
stroys the unforced part of h(x,t). At the same time, 



the concentration field begins to form domains. At later 
times, when L x (t), L y (t) ~ 2iv/k, the domains align with 
the gradient of the forcing term. The growth of the 
domains continues in this direction and is arrested (or 
slowed down considerably) in the direction perpendicu- 
lar to the forcing. The domains are string-like, with kinks 
occurring along lines where a (x,y) is minimized, as ev- 
idenced by Fig. [2] (a)-(d). The growth of L x and L y is 
shown in Fig. 0[ It is not clear whether L y is arrested or 
undergoes slow linear growth and so we do not report its 
growth rate. 

For moderate values of the backreaction strength with 
r ~ 0(1), the height field again assumes a profile aligned 
with the surface tension, while domains of concentration 
now align in a direction perpendicular to the forcing gra- 




FIG. 3: The height field for r = and t = 30000 aligns with 
the applied surface tension. The height field at t = 30000 for 



1 is similar. 



dient. Domain growth continues in the perpendicular 
direction and is arrested in the direction of the driving- 
force gradient. A pattern of string-like domains emerges, 
with domain boundaries forming along lines where both 
a (x, y) and h (x, y, t) are maximized. Eventually, the do- 
main boundaries align perfectly with the surface tension 
maxima, as evidenced in Fig. [2] (e)-(h). 

The control of phase separation by surface shear there- 
fore depends crucially on the backreaction. This result 
is amplified by the existence of a no-rupture condition 
only for the r = case (no backreaction). This con- 
dition relies on the alignment of the height and surface 
tension profiles, which is exact only when the backre- 
action is zero. Then, at late times, the system evolves 
towards equilibrium and is described by the steady state 
V_l • (\h 2 V '±a — ^h 3 X7±p) = 0, which by the alignment 
property reduces to the one-dimensional equation 
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By multiplying both sides of the expression by h, differ- 
entiating and then evaluating the result at #o, a mini- 
mum of both surface tension and height, we obtain the 
condition 
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Since xo is a minimum of height, h" (xq) > 0, which 
prevents h (xq) from being zero. On the other hand, for 
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r and cfq sufficiently large, the alignment of height and 
surface tension profiles is not exact, the one-dimensional 
state is never reached and hence the result in Eq. (|4]) does 
not apply. In that case, simulations show that the film 
ruptures in finite time. 

Given an applied surface tension gradient, we have out- 
lined, by numerical simulations and calculations, three 
possible outcomes for the phase separation, depending on 
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FIG. 4: (Color online) Growth of L x and L y for (a) r = 0, 
where L x grows and L y saturates or undergoes slow growth. 
Since the growth rate of L y is small, it is not clear whether 
saturation or slow linear growth takes place; (b) r = ^, where 
L x saturates and L y grows. 

the backreaction strength r. For r<l, the concentra- 
tion forms string-like domains, aligned with the applied 
force. For r ~ 0(1), the concentration forms domains 
that align perfectly in a direction perpendicular to the 
applied force. For r ^> 1, the forcing causes the film 
to rupture. The inter facial tension or backreaction must 
therefore be chosen carefully in a real fluid to achieve the 



desired outcome. 

In conclusion, we have derived a thin-film model of 
phase separation based on the Navier-Stokes Cahn- 
Hilliard equations, in which the reaction of concentra- 
tion gradients on the flow is important. We have used 
this model to give a qualitative picture of the features 
of phase separation in real thin films, in particular the 
tendency of concentration gradients to promote rupture 
in the film, and to produce peaks and valleys in the free 
surface that mirror the underlying domain morphology. 
We have found that in the presence of a unidirectional si- 
nusoidal variation in surface tension, the strength of the 
backreaction determines the direction in which the do- 
mains align. This result could prove useful in microfab- 
rication applications where control of phase separation is 
required [111 ]. 

Because the lubrication model suppresses vertical vari- 
ations in the concentration field, we are limited to the 
case where the binary fluid components interact identi- 
cally with the boundaries at the substrate and free sur- 
face. However, the model quite generally gives an ac- 
curate description of surface roughening arising from van 
der Waals forces. More detailed models based on this ap- 
proach, involving different boundary conditions that bet- 
ter reflect wetting behaviour 0, Hl[ and a concentration- 
dependent Hamakar coefficient, will capture a wider 
range of thin- film behaviour. 
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